A Protocol For Cooling and Controlling Composite Systems by Local Interactions 
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We discuss an explicit protocol which allows one to externally cool and control a composite system 
by operating on a small subset of it. The scheme permits to transfer arbitrary and unknown quantum 
states from a memory on the network ( "upload access" ) as well as the inverse ( "download access" ) . 
In particular it yields a method for cooling the system. 
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I. INTRODUCTION 

Repetitive applications of the same quantum trans- 
formation have been exploited to achieve noise protec- 
tion [J . cooling, state preparation @, S 3 , and quantum 
state transfer [5|. Motivated by the above results, in 
Ref. we developed a scheme for controlling larger sys- 
tems when control is only assumed to be available on a 
subsystem. Once this is achieved, apart from cooling and 
state preparation, it is also possible to perform arbitrary 
quantum data processing (e.g. measurements, unitary 
rotations). This is similar in spirit to universal quantum 
interfaces of Ref. [3] , but our approach allows us to spec- 
ify explicit protocols and to give lower bounds for fideli- 
ties. These techniques are also related with the "asymp- 
totic completeness" property introduced by Kiimmerer 
and Maassen 0, Q which allows one to control a system 
by coupling it with quantum mediators. 

In the present paper wc review the scheme of Ref. @ 
by showing how arbitrary quantum states can be written 
into (i.e. prepared on) a large system, and read from it, 
by local control only. This implies that arbitrary quan- 
tum operations on the system state can be performed. 
An important specific task is the cooling of the system 
to its ground state. Using some heuristic argument, we 
will provide an estimate of the convergence time of the 
cooling and we will test it with some numerical examples. 
We develop the protocol in several steps. First, we show 
that the system of interest can be actively brought to its 
ground state by replacing its controlled part with fresh 
"cold" qubits from a memory. We then find that cool- 
ing implies that the information about the initial system 
state is transferred into the memory, and design a lin- 
ear map that decodes this information. Since this map 
is generally not unitary, we use the polar decomposition 
to find its best unitary approximation. The fidelity of 
information decoding can then be lower bounded by the 
overlap of the system state with its ground state. Finally 
we design the reverse operation allowing us to transfer in- 
formation from the memory to the system. 

The material is organized as follows. In Sec. |TT] the 
protocol is presented in its general lines. In Sec. IIIII 
we give a detailed derivation of the coding and decod- 
ing transformations and derive bounds for the fidelities. 
Numerical estimations of the protocol performances are 
given in Sec. IIVI focusing on the case of locally controlled 
Heisenberg-like coupled spin networks. Conclusion and 



remarks are in Sec. |V] while technical material is pre- 
sented in the Appendices. 



II. THE PROTOCOL 

Consider a composed system described by the Hilbert 
space Ti ~ Tic ® T^c ® ■ We assume that full control 
(the ability to prepare states and apply unitary transfor- 
mations) is possible on system C and M, but no control 
is available on system C. Moreover, we assume that C 
and C are coupled by a time-independent Hamiltonian 
H. We show here that under certain assumptions, if the 
system CC is initialized in some arbitrary state we can 
transfer ("download") this state into the system M by 
applying some operations between M and C only. Like- 
wise, by initializing the system M in the correct state, 
we can transfer ("upload") arbitrary states on the sys- 
tem CC. The system M functions as a quantum memory 
and must be at least as large as the system CC. As 
sketched in Fig. [T] we can imagine the memory to be split 
into sectors Me, having the same dimension of C, i.e. 
Ti-M — ®e=i T^Me with dimT^Mf = dimT^c- 



A. Downloading info from CC to M 

The downloading protocol we present here is composed 
by two stages: a swapping stage, in which at regular time 
intervals wc couple the subsystem C to the first L memo- 
ries M; and a decoding stage in which wc apply a unitary 
transformation to the first L memories in order to re- 
cover the initial state of CC. As we will see for any finite 
L our analysis does not guarantee that the fidelity be- 
tween the recovered state and the initial state of CC is 
perfect. However, in Sec. IIII A1 it will be shown that by 
augmenting L one can make the fidelity arbitrarily close 
to one. 

We assume that the memory is initialized in a factor- 
ized state of the form 



|0) 



M 



|0) 



M, 



(1) 



where |0) is a state whose properties will be specified 
in the following. To download a generic state, we let the 
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FIG. 1: Example of the model discussed in the text. Here the 
system CC* is formed by 7 spins characterized by some time- 
independent Hamiltonian H (the coupling are represented by 
red lines connecting the spins). The system CC can only 
be controlled by acting on a (small) subsystem C (in this 
case represented by the uppermost spin of the network) . The 
coupling H can - in some cases - mediate the local control on 
C to the full system CC. In our case, system C is controlled 
by performing regular swap operations St between it and a 
2-dim quantum memory Mi. 



system CC to evolve for a while according to its Hamilto- 
nian i7, perform a unitary gate which couples C to one 
of the sectors of Af , let CC evolve again and so forth. 
More specifically, at step i of the protocol we perform a 
unitary swap Sg = Scm^ between system C and system 
Mi 0. After the Lth swap operation the protocol stops. 
This is the swapping stage and it is characterized by the 
unitary operator 

W = SlUSl-iU ■■■SgU ---SiU , (2) 

where U = exp [—iHt\ is the time-evolution of CC for 
some fixed time interval t. As discussed in Ref. @, the 
reduced evolution of the system C under the transfor- 
mation (21) can be expressed in terms of the following 
completely positive CP map 

T{pc) = Trc[U{pc®\Q)cmU^] , (3) 

where lO)^ is the state that is swapped in from the mem- 
ory and Trc[- • • ] indicate the partial trace over the sub- 
system C. Indeed, after L swaps the state of C is ob- 
tained by taking the partial trace with respect to C and 
M of the vector W (|^/')cc|0) m) where \i}))cc the initial 



state of CC, i.e. 

pg^' = TrcM [W^(|^')cc(*I<»|0)m(0|)VF^] 

= TOTo---or (p'g) ^ T^-\p'c) , (4) 
L — I times 

where "o" represents the composition of super- 
operators and 

p^^Trc[C/|V^)cc(V'|t/t] , (5) 

(for an explicit derivation of this expression see Ap- 
pendix \Q . 

Our main assumption is that the map t is ergodic with 
pure fixed point which we denote as \Q)c. Explicitly this 
means that the only state which is left invariant by r is 
the vector |0)p, i.e. 

r{Pc)^Pc ^ Pc-|0)c(0|. (6) 

As shown in Refs. [ll|, this implies that the channel 
T is relaxing (mixing), that is 

lim r"(ac) - |0)c(0| , (7) 

n — >oc 

for all Uq. This condition gives rise to the controllability 
of the system. Indeed from Eq. ([7]) it follows that for 
sufficiently large L, an initial state of the form |'(/')cc ® 
|0)m can be approximated as 

W^(l^>cc ® |0>A/) « |0)cc ® \^mM. (8) 

The right hand side of this equation factorizes into pure 
states because the transformation W is unitary, and be- 
cause both the initial state of CG and M and the final 
state of CC are pure. This implies that, in the asymp- 
totic limit of infinitely many protocol steps (i.e. L 3> 1), 
the system CC has been "cooled" into the state \S)cc 
while all the information regarding the initial state \tl!)cc 
must be contained in the vector |<i>(i/)))M [l^. Further- 
more, it is at least intuitively clear that such information 
can be recovered by the application of a proper unitary 
"decoding" operation V"^ on M which does not depend 
on the input state of the system {decoding stage), i.e. [l3| 

|$(^))a/ « . (9) 

At a mathematical level, the convergence of the down- 
loading protocol described above only depends upon the 
invariant property ^ — see Ref. 0. In Sec. |TTT]we will 
briefly review such a proof and provide a characterization 
of the unitary transformation V. 

B. Uploading info from M to CC 

For uploading states on the system CC, we again make 
use of the unitarity of W. Let us again first give a simple 
hand-waving argument why this is possible. 
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Suppose you want to drive the system into the state 
\'4')cc- To do this, you first use the downloading proto- 
col to make sure that the system is in the state 10)^(5 
("cooling"). Then you bring the memories into the state 
I $(■)/')) 7\/ they would have been ended up in case one was 
trying to download \4')cc from CC into M as in Eq. ([8]). 
Now the quantum recurrence theorem [Toj implies that 
there is a m such that 

I^>cc®|0)a/ « W^{\4,)cc®\^)m) 

« W"^-WQ)cc®M^))m) , (10) 

where we have made use of Eq. ([8]) . Hence by applying 
W m times you have approximately initialized \'4')cc- 
course it remains to be shown that unknown states can 
be written to the system, too. This and the mathemati- 
cal details will be discussed in the next section. Another 
problem with Eq. (|10p is that the recurrence parameter m 
typically needs to be huge^ scaling double exponentially 
with the number of qubits in the system. There are how- 
ever alternative, more efficient ways of implementing an 
uploading process from AI to CC. The simplest one is of 
course to apply the inverse transformation W^^ = 
to the state of Eq. ^ . Indeed the protocol we presented 
in Ref. m is based on this idea, which is a generaliza- 
tion of [4!|. Unfortunately the inverse of W is generally 
unphysical in the sense that it requires backward time 
evolutions U~^, i.e. one would have to wait negative 
time steps between the swaps (see however Ref. [15| for 
cases in which such an inverse time evolution can be im- 
plemented by clever external control techniques). 

To overcome this problem we introduce an extra hy- 
pothesis. Specifically we consider the case in which the 
invariant property ^ holds also for the channel r' ob- 
tained from Eq. ^ by replacing U with W , i.e. 

T'ipc) = Trc [U^ {pc ® |0)c(0|) U] . (11) 

Under this condition, similarly to the case of W discussed 
in the previous section, one can verify that in the limit 
of large i) the transformation 

W = SlU^Sl-iU'' ■ ■ ■ SeU'< ■ ■ ■ SiU^ , (12) 

applied to |V')cc |0)m will converge to a vector of the 
form \0)cc ® I^'('/'))m; i'i) there exists a unitary trans- 
formation V' which does not depend upon lip) and which 
applied to AI gives 

v'^<i>'{i,))M - mM ■ (13) 

From this we can write 

Wcc ® \0)m « {Wy V {\0)cc ® Wm) 
« {USi ■■■USf- USl-iUSl) V {\Q)cc ® I'/'>a/)- 

(14) 

What it is relevant for us is the fact that now the unitary 
transformation on the input state \Q)cc® IV') a/ does not 




FIG. 2: Summary of the scheme: any CP map A can be ap- 
plied to the system by acting on the memory instead through 
the transformations shown in the figure. The red and green 
areas represent the downloading and uploading part of the 
protocol, respectively. The unitary operators W and W'^ of 
Eqs. ((2} and (|15p are generated by acting on the memory and a 
small subsystem of the system only; and V' are instead the 
decoding and encoding unitary transformations introduced in 
Eqs. Q and (|13|l . respectively — see also Sec. IIIII 



involve "time-reversal" evolutions f/"^ but only "proper" 
time evolution U . Therefore, by imposing the condi- 
tion ([6]) on r', we are able to define an uploading pro- 
tocol which transfers an unknown state \'ip) from M to 
CC . Similarly to the downloading scheme it is composed 
by two stages: an encoding stage in which we apply the 
unitary transformation V to "prepare" the memory M 
and a swapping stage in which we apply the unitary 

{W')^ = USi---USf-USL-iUSL . (15) 

by recursively coupling C to the M through swaps. 

Two remarks are mandatory. On one hand, as in the 
case of the downloading protocol, the convergence of the 
transformation (|14p only depends upon the invariant con- 
dition ^ of the channel t'. On the other hand, there 
exists a large class of physically relevant Hamiltonians H 
(e.g. nearest neighbors Heisenberg coupling Hamiltoni- 
ans) for which both r and r' verify the such condition — 
we refer the reader to Ref. 01 for details. For such Hamil- 
tonians, our analysis will yield both a simple downloading 
and uploading mechanism. Putting these two elements 
together one can also realize more sophisticated controls. 
For instance, as shown in Fig. [21 one can perform any 
quantum transformation A on CC by first downloading 
its state on M, transforming it, and finally uploading the 
final state back into the system. 



III. CODING TRANSFORMATION 

In this section we derive the decoding transformation 
yl' that relates states on the memories Ad to the states on 
CC in the downloading protocol. To do so we exploit the 
formal decomposition of the evolved state of the system 
after L steps (see Appendix |B|) . The encoding transfor- 
mation V' of the uploading protocol can be obtained in 
a similar way. 

Consider an orthonormal basis {|'0fe)cc'} of '^cc- 
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cording to Eq. (jBl[) after L swaps it becomes 

W{\i:k)cc\^)M) (16) 



|0)c® V^|0)c|(/)fc)M + - Vk\^k)cM 



where |Afc)i5^j is a vector orthogonal to and rjk ~ 1 
as in Eq. (|B6p . This equation shows that with high prob- 
abihty, the transformation W maps the orthonormal vec- 
tors \ipk)cc i'^lo ^^^^ vectors \(t>k)M of the first L mem- 
ories. For any finite choice of L, the latter are typically 
not mutually orthogonal. However one can use Eq. (|B6p 
to show that in the limit of large L the vectors \4>k)M 
become approximately orthogonal. Indeed from the uni- 
tarity of W and from Eq. ITBl) and (|B3p we can establish 
the following identity 

Ali^klMM (17) 
= Vk' Skk' + \/vk (1 - Vk') CCM{'4'kO\Ak')cCM 

+ VVk' (1 - Vk) i5CM{^k\^Jk'0)cCM 
+ V (1 - ??fc)(l - Vk') cCM{^k\Ak')cCM ■ 

To simplify this expression we define rjo = miufe rjk ■ Since 
Eq. (jB6p applies to all r]k the parameter r/o must satisfy 
the inequality 



l-rio^^K {L- ifo K 



L-l 



(18) 



Furthermore from Eq. (fTT]) it follows that for k ^ k' one 
has 



A/ 



(</>fc|(/>fc' 



WVk (1 - %') lcCA/(V'fcO|Afc')ccMl 
+ v/f?fc' (l-77fe) lcCA./(Afe|V'fc'0)(5cA/l 
+ V(l-'7fe)(l-^fe') lcCA/(Afe|Afc')cCA//l 

s$ 2v/l - ?7o + (1 - ??o) 3^1 - Vo , 



(19) 



which according to Eq. and using the fact that the 
parameter k is strictly smaller than 1, shows that for 
large L the vectors \(t>k)M and \4>k') m become orthogonal. 

Define then the linear operator D on 7ij\/ which per- 
forms the following transformation 



D\lpk)M = |'?!>fe)A/ 



(20) 



with \ijjk)M being orthonormal vectors of M which rep- 
resent the states {|'0fc)cc} of ^cc- Formally they are 
obtained by a partial isometry from CC to M and are 
"good" representations of the \'4'k)cc- The operator D 
in some sense "corrects" the non-orthogonality of the 
|0fc)A/: indeed its inverse (when definable) allows us to 
pass from these approximate images of the \'ipk)cc lo the 
good representations \ipk)M- Therefore seems to 

be a good candidate for defining our decoding transfor- 
mation V. Unfortunately however D is NOT unitary (it 
maps an orthonormal set of states into a non-orthonormal 
one) and typically will not be even invertible. 



The idea is then to replace D with its best unitary 
approximation V [itI . p 432]. The latter is obtained by 
taking the polar decomposition of D, i.e. 



D = PV , 



(21) 



with P positive semidefinite. The unitary V minimizes 
the norm distance from D yielding the inequality 



\D^V\\2 




^ V3 dec (1 ' VoY^^ 



(22) 



where we introduced the eigenvalues Afe of D'' D and 
used Eq. to bound them according to the inequal- 
ity I '^fc ^ 1| ^ 3 d(j(j y/1 — 770 (in these expressions 
d(j(j = dcd(j is the dimension of the system CC and 
||e||2 stands for ^/J2kk' IQfcfc'P with Qkk' being the ma- 
trix elements of the operator 8). The inequalities ([22]) 



and (US]) show that for L 00, D can be approxi- 
mated arbitrarily well by the unitary V. We can hence 
define V'' as our decoding transformation which inverts 
the mapping (|20p and transforms the "bad" representa- 
tions |0fc)Af of the \'ipk)cc i'^lo "good" representa- 
tions {ipkjM- It is worth stressing that, by construction, 
V does not depend upon the input state \'ip)cc the 
system CC. 

As mentioned in the introduction of this section, a sim- 
ilar procedure can used to defined the encoding protocol 
of the uploading protocol. Without entering into the de- 
tails we simply notice that in this case D and the vectors 
\(l)k)M will be defined by replacing W of Eq. (fTB)) with 
the transformation W' of Eq. (fT2)) . Taking the polar de- 
composition of such new D it will yield the unitary V' 
which will be used as encoding for the uploading scheme. 

In the following section we will evaluate the transfer 
fidelities associated with such a choice of decoding and 
encoding transformation, showing that they can arbitrar- 
ily increased by choosing L sufficiently high. 



A. Fidelity of the downloading protocol 

Let |V')cc ~ J2k'^k\'>Pk)cc be a generic input state 
of CC. To evaluate the downloading fidelity Fdown as- 
sociated with our decoding scheme we need to compare 
the state of M at the end of the protocol with the state 
\^)m = J2k (^k\ipk)M, i-e. 



own 



(23) 



Here V'' is the decoding transformation defined in the 
previous section, and Rm is the state of the memory 
after the application of the unitary W, i.e. 

Rm = Trcc[M/(|V)cc(V^I®|0)A/(0|)Tyt] 

= ^|0)a/(0| + (1-?7)'^a/. (24) 
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In the above expression we used Eqs. (jBip and (|B2p and 
introduced the density matrix um = Tr(j[| A)(;jj,f (A|]. By 
hnearity we get 

^down(V') = ^ \M{cp\V\i,)M? 

> 7? \M{mi^)M? ■ (25) 

We now bound the term on the right hand side as follows 

|m(0|v^IV')mI = |a./(0I^-^ + ^IV'>a/I (26) 

^ \M{cl)m)M\~\M{(b\D ~V\ij)M\ , 

and use the inequality ([^ to write 

\m{<P\d-v\^p)m\ ^ \\D-v\\2^ Video (1-%)'/* • 

If \i')M was a basis state \tpk)M, then a/ | = 1 

by the definition Eq. (|20p of D. For generic instead 
we can use the linearity to find after some algebra that 

\M{(l>\D\i')M\ >Vm - ^dcc V^-Vo ■ (27) 
Therefore Eq. (|26p gives 

which replaced in Eq. ((25|) yields 

^^dow„(V') > m - 10 rfcc (1 - VoY^^ - (29) 

for all input states \ip)cc- According to Eq. p8)) it then 
follows that by choosing L sufficiently big our download- 
ing protocol will yield transferring fidelities arbitrarily 
close to one. 



B. Fidelity of the uploading protocol 

Following the analysis of Sec. (jll Bp the fidelity for up- 
loading a state \tp)M into CC is given by 

F,p(V.) = (30) 
cc(V'|TrA/ [W'^V (|V)a/(V'I ® |0)cc(0|) V'^W] H')cc- 

A lower bound for this quantity is obtained by replacing 
the trace over M with the expectation value on |0)a/, i.e. 

FupW ^ cc(V'|a/(0|M^'V(|V')a/(^|® |0)oc(0|) 

cc 

= \cc{0\M{m'^W'\0)M\i^)cc\^ (31) 

In deriving this equation we used Eq. (jlSp and a decom- 
position of the form of Eq. (|Bip to simplify the vector 
W^'|0)a^ |V')cc- this case ?/ is defined as in Eq. (|B5p 
with r being replaced by r' of Eq. pip . Since we are 
assuming that this CP map satisfies the condition ^ 



it follows that also rj' obeys an inequality of the form 
(|B6p with K and k replaced by new constants K' and 
k' g]0, 1[. We also notice that last term of Eq. ([?T|) has 
the same form of the lower bound (P5)l of the download- 
ing fidelity. Therefore, by applying the same derivation 
of the previous section we can write 

F^pii;) ^ rj'o - W dec a - Vo)'^\ (32) 

with 

This shows that, as in the downloading case, also the up- 
loading fidelity converges to unity in the limit of large L. 

IV. EFFICIENCY OF COOLING 

In this section we provide some numerical estimation 
of the quantities 77 of Eq. (|B5p which measure the prob- 
ability of finding the state C in lO)^;. As seen in the 
previous sections this is the fundamental parameter to 
bound the fidelities of both the downloading and upload- 
ing protocol. Moreover, given an initial state \tp)ccj V 
measures the success probability of "cooling" it down to 
the state \0)cc during the downloading process. Accord- 
ing to Eq. (|B6p the quantity rj will asymptotically con- 
verge exponentially fast to unity. However Eq. (|B6p does 
not tell us from which point onwards the convergence is 
exponentially fast, so it would be nice to have alternative 
ways to estimate the convergence speed. 

To simplify the analysis in the following, we will con- 
centrate on the spin network model of Fig. [1] assuming a 
Heisenberg Hamiltonian of the form 

H= J2 d,,AX,X,,+Y,Yy+Z,Zy), (34) 

(ij')6G 

which conserves the total magnetization along the z axis 
(here Xj, Yj and Zj are the Pauli operators of the j- 
th spin and the summation is performed over all the 
edges of the weighted graph G associated with the net- 
work). Moreover we will take the vector |0)c to be the 
configuration where all the qubits of C are in the spin- 
down state, i.e. |0)c = |00---00)c- For this choice of 
the controller state our main assumption of ergodicity 
Eq. ^ is numerically found to be correct for the cou- 
pling graph depicted in Fig. [T] The fixed point is given 
by \0)(j = |00 • • • 00) c (more general conditions of ergod- 
icity for Heisenberg models are given in Q). In this con- 
text 77 coincides then with the probability Pg^'' of finding 
no excitations on the system after L steps of the proto- 
col. Some numerical examples showing the dependence 
of 77 upon the initial state are presented in Fig. [31 As 
expected, asymptotically Pg^'' is seen to converge expo- 
nentially fast. 

An approximate estimation of Pg^'* can be easily ob- 
tained by looking at the average number of spin-up on 
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Number o1 swaps 

FIG. 3: Convergence of the cooling protocol for the weighted 
graph of Fig. [l] where the couplings among the spins is 
given by the Hamiltonian (|34|) (the values of constants dj j/ 
have been chosen to be proportional to the length of the 
graph edge). Four different states \'4')cc ^'i'^ considered: the 
fully flipped state |lllllll)cc, a GHZ state (| 0000000) cc + 
|lllllll)(j(;j)/\/2, a fully mixed state Pqq — 1q(j/2^, and a 
W state ^ (1 1000000) cc + 1 0100000) cc5 + • • • + 1 0000001) cc). 




5 10 15 20 25 30 35 



FIG. 4: Comparison of the approximation Eq. (|37|) with ex- 
act numerical results. Shown is the average number of excita- 
tions (N) on an open Heisenberg spin chain with 7 sites 
\ / cc 

and equal couplings as a function of the number of swaps to 
the memory. The initial state is taken to be | lllllll)(;j(j , 
i.e. with a maximal number of excitations. The three curves 
correspond to different sizes of the region \C\ controlled by 
Alice, and the time interval t has been chosen for each curve 
independently to fit the approximation given in Eq. (|37p . 



CC after L swaps, i.e. 



N 



Tr 



CC 



cc 



(35) 



with p^^l = TrM[VF(|V)cc(V'l ® |0)M(0|)VFt] being 
the reduced density matrix of CC and with N = 
X^fcec c (^'s "t" ^) (here Zk is the z-Pauh matrix of the 
A:-th spin). These quantities are related by 



Pr 



> 1 



N 



(L) 



CC 



(36) 



To get an approximation for the average number of exci- 
tations on the graph we assume now that the time inter- 
val t is chosen such that U shuffles the excitations on the 
graph in a fully random way. For specific systems and 
specific times intervals, this "classical" behavior might 
not be true due to interferences, but for general times 
it is a good approximation (sec Fig. [4]). Let |C| be the 
number of edges on the graph controlled by Alice, and 
\C\ the number of uncontrolled edges. On average, each 
swap takes approximately a ratio |C|/(|C| -I- \C\) of exci- 
tations from the graph to the memory. We then get 



N 



CC 



N 



(0) 



1 



cc \l + \C\/\C\ 



(37) 



This is a reasonable result which shows that the fidelity 
depends on the initial number of excitations and on the 
relative size of the controlled region with respect to the 
uncontrolled region. 



V. CONCLUSION 

We have given an explicit protocol for controlling and 
cooling a large permanently coupled system by accessing 
a small subsystem only. As we have shown, the applica- 
bility relies only on the invariant property ([6]) of a CPT 
map. Since we had to assume a large quantum memory 
in order to control the system, this protocol is not useful 
for replacing control in a homogeneous setup, but may 



coherent 
lightsource 




^ ^ 



readout to QC 

FIG. 5: A CCD-like application of our protocol could allow a 
light sensitive array of qubits to be read out coherently by a 
quantum computer without "disturbing" the qubits much. 
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well have applications in inhomogeneous scenarios (when 
control is harmful or expensive in some regions but easy 
in others). For example, we imagine a CCD-like appli- 
cation, in which a set of permanently coupled qubits is 
read out by a Quantum Computer in a coherent manner 
(see Fig. [5]). 
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Here we derive the evolution ([4]) of the subset C in 
terms of the CP map ([3]) . First rewrite the reduced den- 
sity matrix ([4]) as follows 



pjf^ = TrcM[W^(|«')cc(*l®|0>M(0|)W^t] 



= Tr, 



c 



• • Tr 



Ah 



Ml 



cc 



Ml 



(Al) 



For the sake of clarity it is useful to explicitly denote the subsystems on which the various operators are acting on 
(e.g. Qab indicates that the operator 6 acts non trivially only on the subsystems A and B, while it is the identity 
elsewhere). By doing so and by using the properties Q of the swap it is easy to verify the following identities: 



Tr,. 



M, 



SiU(^Pc ® |0)c(0| ® \0)mAO\)u^sI 



Tr 



M, 



s, 



'cM,Ucc{pc ® |0)c(0| ® \0)MAO\)ulcShM, 



Tr 



Mi 



ScMiUcc(sljM^ScMi) (pc ® |0)c(0| ® \0)m, 



CMe •^CMe I U^,^,S(j^j^ 



Tr 



Mi 



UMic{pc ® |0)Af,(0| ® |0)c(0|)t/],^J = t{pc) ® |0)c(0| 



(A2) 



which holds for all pQ and Equation ^ then follows by replacing this into Eq. (|Aip for all ^ > 2 and by using the 
identity 



TrA/^ 
= TrA/; 



SiU{\^)ccm®\Q)Mi{Q\)u^S\ 



Tr 



Ml 



s, 



MiC 



= Tr 



'cMiUcc{\^)cc 

®\0)c{0\=p'c 



(A3) 



U 



MiC 



MiC 



MiC 



with p'^ as in Eq. (O. 



APPENDIX B: DECOMPOSITION EQUATIONS 

Here we give a decomposition of the state after ap- 
plying the W operator of Eq. ([2]) . This will allow us to 
estimate the fidelities for state transfer in terms of the 
relaxing properties of the map t. 

Let \'ip}cc € T^cc be an arbitrary state. We notice 
that the C component of M^|V')cc|0) a/ is always |0)c- 
Therefore we can decompose it as follows 



W\^P)cc\0)m = |0)c ® [V^|O)c|0)a/ + V^\^)cM 

(Bl) 

with rj G [0, 1] and with | A)(jAf being a normalised vector 
of C and M which satisfies the identity 



c(0|A)cA/ = 



(B2) 



It is worth stressing that in the above expression t], |(/))a/ 
and I A) (J A/ are depending on \ip)cc- ^ similar way we 
can decompose the vector obtained by acting with W"^ 
on the first term of Eq. (jBl|) . i.e. 



W^\0)cc\^)m = |V')cc|0)m + |A)ccM, 

(B3) 

where |A)p(5A/ is the orthogonal complement of 



cc 



(V'|m(0|A) 



CCM 



= 



(B4) 



Multiplying Eq. p33)) from the left with cc (V'|Af(0| and 
using the conjugate of Eq. (|B1[) wc find that f] = ij. An 
expression of 77 in terms of t can be obtained by using 
Eq. Therefore from Eq. (|B1[) and the orthogonality 
relation (IB2D it follows that 



V = c{ 



-L-l 



p'c) I0)c, 



(B5) 



8 



which, since t is rela xing: , shows that 77 1 for L 00. where = ■\/Tr[8t0] is the trace norm of the opera- 

Moreover we can use [3| to claim that tor 0, K is a constant which depends upon dc = dimTig, 

and where k e]0, 1[ is the second largest of the moduli of 
l-V = Ic(0|t^"^ {p'c) |0>c - 1| eigenvalues of r. 

< ||r^-^(pb)-|0>c(0|||i 
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